The Mass of the Cosmos 
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We point out that the mass of the cosmos on gigaparsec scales can be measured, owing to the 
unique geometric role of the maximum in the areal radius. Unlike all other points on the past null 
cone, this maximum has an associated mass, which can be calculated with very few assumptions 
about the cosmological model, providing a measurable characteristic of our cosmos. In combina- 
tion with luminosities and source counts, it gives the bulk mass to light ratio. The maximum is 
particularly sensitive to the values of the bulk cosmological parameters. In addition, it provides a 
key reference point in attempts to connect cosmic geometry with observations. We recommend the 
determination of the distance and redshift of this maximum be explicitly included in the scientific 
goals of the next generation of reshift surveys. The maximum in the redshift space density provides 
a secondary large scale characteristic of the cosmos. 

PACS numbers: 98.80.-k, 98.65.-r 



I. INTRODUCTION 



The emergence of automated large-scale redshift sur- 
veys 0, is opening up new possibilities for measuring 
the content and dynamics of the cosmos on very large 
scales. We point out a significant characterisation of the 
cosmos on gigaparsec scales that will become measurable 
with the next generation of surveys. 

It is well known in observational cosmology that our 
past null cone (pnc) has a maximum in its areal radius, 
Rrm where the angular size of sources of a given size 
is minimum. Beyond this point, more distant images, 
though dimmer, subtend larger angular sizes It 
is also known in relativistic cosmology that this maxi- 
mum occurs where the observer's pnc crosses the appar- 
ent horizon. What hasn't been realised is that a measure- 
ment of this maximum is equivalent to a measurement of 
the mass within a sphere of areal radius Rm, and that 
this relationship is quite general, not requiring the as- 
sumption of homogeneity for example. Less well known 
is that the number density of sources in redshift space 
also has a maximum. This latter maximum, however, 
does not have such a deep significance as the former. 

It is best to use the Lemaitre-Tolman (LT) model, as 
both R and M are primary functions, so the relation- 
ships we are interested in are particularly clear. Also, 
the equation of state — pressure free matter plus A — 
is very suitable for the post-recombination universe, and 
the spherical symmetry of the model about the origin 
is entirely natural in the context of observations, since 
isotropy is fairly well established, and of course our own 
past light cone is centered on ourselves. 



The Lemaitre-Tolman (LT) metric 51, is 



ds^ = -de 
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where R ~ R{t,r) is the areal radius, R' — dR/dr, and 
E = E{r) is an arbitrary function of coordinate radius 
r. We use geometric units in all equations. This metric 
describes pressure-free matter in comoving coordinates. 
From the Einstein field equations. 
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where k = SttG/c* = Sir, p{t,r) is the density, R — 
dR/dt and M = M(r) is also arbitrary. The solutions to 
are best found numerically for the general case. The 
density is divergent at the bang and crunch, i? = 0, and 
also at shell crossings, where i?' = but M' ^ 0. The 
latter singularity is avoidable but the former isn't. 

As can be seen from , the function E determines the 
local geometry. In addition, (0) shows that it is a kind of 
energy parameter. Function M is the gravitational mass 
within a comoving sphere of radius r, so it includes any 
putative dark matter component (s), but does not include 
the "density" associated with A. This M is not the same 
as the integrated proper density M, 
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pd^V^ / (M'/x/l + 2E)dr 
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This latter is the mass one would obtain by summing the 
masses of individual galaxies, gas clouds, dark matter 
concentrations, etc. See for further discussion of the 
LT model and supporting references. 
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II. OBSERVABLES 

The light rays making up the pnc (past nuU cone) are 
the incoming radial null geodesies arriving at the central 
observer at a given moment, in particular: 
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At 



Ar ^\ + 2E ' 



t = to at r = 



(5) 



with solution t = t{r). Any quantity Q(t,r) evalu- 
ated on our current pnc will be indicated with a hat: 
Q — Q('") — Q{iiT)Tf)i a-iid for expressions a square 
bracket with subscript "n" will be used. The redshift 
of an observed object, located at rg, is given by 
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where R' may be found from l^l- Along a ray, the areal 
radius is i? = R{t{r),r), and its rate of change along the 
ray, using Q and (O, is 
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Since the coordinate r is not observable or physically 
meaningful we calculate 
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0, so 



Near the origin we have M 0, E ^ and z 

AR/Az^[R'/R%o ■ (9) 

See 0, 01 for further details. 

If a source such as a galaxy has measured angular di- 
ameter S and known actual diameter D, then the diame- 
ter distance is identically the areal radius on the pnc and 
is defined by 



R 



D/S 



(10) 



If the apparent luminosity of an observed source is 
and its absolute luminosity is known to be L, then the 
luminosity distance is defined by 



dr. = 



(11) 



where da = 10 pc, m is the apparent magnitude and ifi 
the absolute magnitude. For a general curved spacetime 
and arbitrary motion, the reciprocity theorem tells us 



dL = doil + zf 



(12) 



li do = R is calculated from m k, z measurements using 
l|12|l . the fractional error is comparable with that of d^, 
assuming the redshifts are fairly accurate: 



5R/R = 0.2 In 10 5{m - m) - 25z/(l + z) 



(13) 



For each of these distances, it is essential to know an 
intrinsic source property, D or L (or m). While their 
values for nearby sources can be determined using other 
distance measures, for distant sources there is the ex- 
tra difficulty of determining how much they have evolved 
with time, or even what type of nearby object the source 
corresponds to. 

If n is the number density of sources in redshift 
space(/steradian/unit redshift interval), and /i is the 
mass per source, then the relation between n and p is 
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See [Tl| for a discussion of multiple source types and mul- 
ticolour observations. 



III. CHARACTERISTIC COSMIC MASS AND 
DENSITY 

In an expanding, non-inflating universe, the diameter 
distance R{r) necessarily has a maximum. The locus 
of such points, where different rays are momentarily at 
constant R, is the apparent horizon (AH). We find it by 
putting i?' = in (T)), giving 



ARl - 3Rra + 6A/„ 



0, 



(15) 



so thus if Rm is on the AH, then 2Mm — Rm — Ai?^/3, 
and if A = this is the familiar R — 2M. In general the 
locus of the pnc, and the variation of measurables down 
it, are strongly affected by the details of the cosmological 
model. This point, however, where the pnc crosses the 
AH, has a unique meaning: 

(a) there is a simple direct relationship between R and 
M that is independent of any inhomogeneities in 
E, tB and M, 

(b) since R is measurable, the mass M within that ra- 
dius is immediately determined, 

(c) this point flags a major causal feature of model, 

(d) the maximum in i? is a distinctive feature of the 
R-z plot. 

No other point on the pnc has such a simple, generic 
M-R relationship, and this holds true whether not the 
model is homogeneous. Therefore Rm and the redshift 
Zm where it occurs are distinctive characteristics of the 
cosmological model. 
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For Mm > & Rm > 0, equation (|15|l has solutions if 
M,„ < 1/(3\/A) = (Mm) max 1 

(16) 

and this maximum possible value of Mm occurs at 

(-Rm)ma. - 1/\/A . (17) 



Thus if an incoming ray reaches R > {Rm)max, then it 
has no maximum in R. Worldlincs with M > {Mm)max 
never meet the AH, and those with M < {Mm)Tnax cross 
it twice. 

Condition (|15|l will always hold somewhere for any cos- 
mology with a non-zero matter density, because Af = at 
the observer, and increases outwards. In all LT models 
with a bang, the AH goes as i? « 2M near the bang. 
In ever-expanding models, the AH asypmtotically ap- 
proaches the de Sitter horizon, R = -\/3/A , as t — )■ oo. 
The "incoming" ray that is tangential to i? = -^/S/A 
at t = oo divides rays that have a maximum and reach 
the origin, from those that have no maximum and have 
ever-increasing R. In closed re-collapsing models, putting 
R = and _E = — 1/2 in ||2J) reproduces (|15|l . thus show- 
ing that the moment of maximum expansion at the max- 
imum in the spatial sections lies on the AH. (In fact it is 
where the past and future AHs cross.) Thus all pncs in 
(physically well-behaved) closed re-collapsing LT models 
have maxima in their areal radii. For further discussion 
of the apparent horizon, see 0, . 

It is evident from H14() that the redshift space density 
also has a maximum, because, although M' > and 
increases from zero at the origin both R' and (1 -I- z) 
are finite at the observer, and diverge towards the bang. 
For the maximum in fin, we put 
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(18) 



and solve for /In„j and Zm- The locus of this maximum 
doesn't have a deep geometric or physical meaning, but 
depends on the redshift behaviour down the pnc, which 
depends strongly on the details of the model. Neverthe- 
less, one can say that the maximum in fin is a large scale 
characteristic density of the model. 



IV. ROBERTSON- WALKER CASE 

The key results are those given above in section IIIII 
which apply to a fairly general class of realistic post- 
recombination cosmologies. However, it is useful to look 



at the characteristic cosmic mass and densitiy in the stan- 
dard homogeneous model. The above results are spe- 
cialised to the homogeneous case in appendix 

The Robertson- Walker (RW) version of the AH equa- 
tion (jA12|l seems quite complex, and hides the key rela- 
tionship that is evident in H15|l . 

The dependence of Rm, Mm, ■M.m and Jm. on Ho, f^A 
and Vim is shown in figQ and for reference the R{z) and 
fin{z) curves for a range of values of these parameters are 
given in the appendix^. 




5- n =0.5 



= 

+10"=70;Omm = 0.3,"OmL = 0.?- " 




^ Here we use a well behaved r coordinate, such as one with M oc 
r'' or R(t = const, r) oc r. 



^ For a related earlier treatment with A = 0, see |14| . 
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FIG. 1: (a) The dependence of Rm and Zm on the RW pa- 
rameters; (b) the dependence of Mm (solid hnes) and A4m 
(dashed hnes) on the RW parameters, (c) the dependence of 
/m^ and z,„ on the RW parameters. 



V. DISCUSSION 

The properties outlined in section 11111 below eq H15|) 
make the maximum in i? a very significant feature of a 
cosmological model, both theoretically and observation- 
ally. 

Both Rm and Zm may be measured, and calculating 
M,„ from H15|) requires only the value of A, and is inde- 
pendent of the details of the cosmology, or even whether 
large-scale homogeneity exists. Note also that Mm is not 
very sensitive to uncertainties in f^A ; for a given measured 
Rm, an increase from to 0.7 decreases the calculated 
Mm by 12%3. 

Current galaxy redshift surveys only extend to z = 0.3, 
while quasar redshift surveys extend to z = 6. However 
the data is not complete enough, and the quasar popu- 
lation too diverse to be useful. The next generation of 
galaxy surveys will extend past the redshift of the maxi- 
mum, and recent supernova observations already 
approaching it. The R-z plot may be obtained from di- 
rect measurements of redshifts and angular diameters, 
or derived from m-z measurements. While the former 
would be ideal, the latter is acceptable, since the reci- 
procity theorem is very general. Of course, knowledge 
(or assumptions) about true diameters and absolute lu- 
minosities of sources, and their z evolution, is essential, 
and at large z this is a significant uncertainty. 

While locating Rm reliably requires a good sample of 
data points in a range near Zm, determination of n{z) is 
not so easy, since one must be sure of detecting or reliably 
estimating all masses. Nevertheless, a sufficiently deep, 
complete survey should provide an indication of /m^ and 

• 

An issue for future consideration is the effect of inho- 
mogeneities on the uncertainty in Rm and Zm determina- 
tions. For example ^3 showed inhomogeneity can cause 
loops in the R-z curve near the maximum. 

The mass-to-light ratio is difficult to determine outside 
gravitationally bound systems, and direct estimates of 
(gravitational) mass from galaxy and cluster dynamics, 
and from gravitational lensing only extend up to cluster 
or supercluster size, whereas the determination of Mm 
reaches several orders of magnitude larger. 

The possibility of determining the metric of the cosmos 
from observations of redshifts, luminosities (or angular 



^ This is the uncertainty in Mm once Rm has been measured, 
rather than the variation in the Mm predicted by a variety of 
models with different values, in which Rm also varies. 



diameters) and the number density of sources, combined 
with the evolution of absolute luminosities, true diame- 
ters, and mass per source, was considered in I^ Hllll^ . A 
project to begin implementing this is now underway |l7l| . 
The cosmic mass Af,„ provides an important cross-check 
on the summed mass at that radius. In fact the theorem 
of — that, given any reasonable set of observations 
and any reasonable source evolutions functions, an LT 
model can be found that fits them — needs to be qual- 
ified, as the combination of observations and evolution 
functions must mesh correctly close to Rm- 

Within the dust-A-RW model, it is noteworthy that 
the maximum is where the R{z) curve is most sensitive 
to variations in and Hq, and nearly so for ri,„, and 
that variations in these 3 parameters move the (z„j, Rm) 
and (primj ^m), loci in very different directions. Thus a 
determination of Rm & Zm would provide rather generic 
limits on Hq, flm and JIa, and combined with the initial 
slope of the R-z graph, or with measurements of /m„ 
and Zm, would fix all three values. Determining Rm and 
Zm to within 10% would by itself provide confirmation 
(or otherwise) of current parameter estimations, while 
5% accuracy would put new constraints on the possible 
values. 

The problem of averaging in GR means that identify- 
ing the RW model that best fits the observations is not 
a well-defined execise. Therefore measurments of bulk 
effects are particularly important, and since Rm, unlike 
any other point on the pnc, gives the total mass on that 
scale, the Rm and Zm values provide a natural definition 
of the best fit RW model. 

Although particular models always have a R-M rela- 
tion, such as IjAip & ljA2p would give for RW models, this 
relation is very model dependent, whereas Mm is not. It 
is more general even than the LT model, as any cosmol- 
ogy will have a locus where the past null cone crosses the 
apparent horizon, and an associated mass is naturally 
defined there. 



VI. CONCLUSION 

In summary, the maximum in the diameter distance 
is the only point on the past null cone that corresponds 
to a model independent mass, thus allowing direct mea- 
surement of a characteristic cosmic mass on gigaparsec 
scales. Therefore it also provides a very large-scale check 
on the mass to light ratio, as well as a reference point 
for determining geometry from observations. Since it is 
a point on the apparent horizon, a measurement of this 
maximum may actually be the first detection of a rela- 
tivistic horizon. 

For RW models, the region near the maximum is where 
R{z) is most sensitive to the values of the RW parame- 
ters. 

We advocate that, with the next generation of surveys, 
direct measurements of angular sizes and hence diame- 
ter distances be compiled and calculated independently 
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of luminosity distances, and that, apart from fitting a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) model 
to the available data, a separate determination of Rm 
be done with a limited data set near thus giving a 
model independent value for M„i. 
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FIG. 2: The R{z) curves for a range of RW parameters: 
(a) different Q,m values {Hq = 70, Q,k = 0.7), (b) values 
[Ho = 70, f7„ = 0.3), (c) Ho values (n„ = 0.3, = 0.7). 
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FIG. 3: The ^,n{z) curves for a range of RW parameters (note 
that the z scale is different from that of fig. @: (a) various 
Q,^ values {Ho = 70, ^a = 0.7), (b) Oa values {Ho = 70, 
Q.m = 0.3), (c) Ho values {Qm = 0.3, Q.a = 0.7). 



APPENDIX A: SPECIALISING THE LT 
EQUATIONS TO RW 

The dust-A-RW model, in its most common coordinate 
system, is obtained if we put 



M = Mor^ , 2E = -kr^ , ts = , (Al) 
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from which we have 

R = rS{t) , = 2Mo/S -k + AS^/3 , 



(A2) 



Ho = So/ So , Kp = (SMo/S^ , Kpo = (iMo/Sl , 

(A3) 

where S is the scale factor. Given Ho^ and ri„i, then 
to = 2/3Ho , npQ = 3n^H^ , A = iH^^A , 



1 — il-A — 



So 

Mo = n,nH^S'y2 



-sign(rife) , 



arbitrary if k = , 

t2 r>3 



(A4) 
(A5) 

(A6) 

(A7) 



where n„i is the baryonic plus dark matter fraction. The 
integrated density equation is 

M=^-^ /sin- Vfc r) ^^Y^ \ (A8) 



so for nearly flat models, where r << 1, 

{M - M)/M « O.Skr^ 

The pnc and AH equations 10), ©, ©, 113 
using Sm = S{t{rm)), become 



VI - kr 



(A9) 
(AlO) 



di?, 



— kr^ — rS 



H{l + z) 



dR„ 



dz 



1 

(All) 

ArlSl - 3r™^„ + 6Afor^ = , (A12) 
6^^o /:; — A , 2/o^ A2n 



kS'3(1 



2rVl - fcr2 S + r^{SS - S^ 





(A13) 



Figs El and 13 show R{z) and p,n(z) for a range of Ho, 
VIa and f2m values. The primary interest is on how these 
parameters affect the maxima. 
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